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I.  INTRODUCTION 


The  finite  element  method  has  been  used  in  the  modeling  of  many  problems 
concerning  elastic-plastic  impact  behavior.  The  correct  geometric  represen¬ 
tation  of  a  mesh  is  important  in  obtaining  accurate  results  without  introduc¬ 
ing  numerical  instabilities  into  the  calculations.  Here,  triangular  elements 
formulated  for  wave  propagation  and  shock  problems  have  been  applied  to  both  a 
classical  elastic  wave  propagation  problem  and  to  the  analysis  of  a  system 
subjected  to  high  velocity  impact  loading. 

The  objective  of  this  study  is  to  ascertain  the  effects  of  different 
orientations  of  the  triangular  elements  and  to  determine  the  optimal  mesh 
size  and  shape  using  EPIC-2,  a  computer  code  for  elastic-plastic  impact 
calculations  in  two  dimensions,  to  best  fit  the  solution.  The  information 
presented  here  is  intended  to  be  of  assistance  to  potential  users  faced  with 
applying  the  finite  element  method  to  impact  problems  and  considering  how  to 
obtain  meaningful  solutions  in  two  dimensions. 

Wave  propagation  in  the  impact  of  a  semi-infinite  elastic  bar  was  chosen 
because  of  the  availability  of  a  closed  form  solution.  The  analytical 

solution  obtained  by  Richard  Skalak,3  and  the  elementary  one-dimensional 

p 

theory  for  elastic  wave  propagation  in  Iona,  slender  rods  were  used  to  Judge 
the  convergence  of  the  numerical  solution  in  the  various  calculations.  After 
obtaining  the  optimal  mesh  site,  the  effects  upon  the  changes  in  artificial 
viscosity  factors  were  observed. 

A  large  deformation  problem  was  also  considered,  in  tne  form  of  a  sphere 
impacting  a  flat  plate  at  high  velocity,  to  asses*  effects  of  grid  size  and 
orientation  on  deformations,  pressures,  and  strains  in  the  presence  of  large 
plastic  flow. 


II.  APPROACH 


Two  numerical  examples  are  provided  to  demonstr  ,e  the  effects  of  mesh 
orientation  and  artificial  viscosity  upon  the  numerics*  solution. 

All  computational  results  were  obtained  with  the  EPIC-2  (Elastic-Plastic 

3 

Impact  Calculation  in  Two  Dimensions)  computer  code,  a  Lagrangian  finite 


1  Skalak,  R. ,  "Longitudinal  Impact  of  a  Semi- Infinite  Circular*  Elastic  Bar," 
Trans.  AStiE,  J.  Applied  Meek. ,  ’/ol  24,  Dec  20 Lb,  pp.  83-63. 

** 

~  Johnson,  W. ,  Impact  Strength  of  Materials ,  Crane,  Runsak  and  Co.,  Hew  York, 
lu7:l,  pp.  J-J2. 

'Jj  Johnson,  U.  R. ,  "EPIC-2,  A  Computer  Prvgi'am  for  Elastic-Plastic  Impact  Cal¬ 
culations  in  Two  Dimensions  Plus  Spin,"  Honeywell,  Inc.,  Defense  Systems 
Division,  Contract  Report  ARBRL-CR-00373,  June  1378  (ADH  A0S8786). 
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element  formulation  where  constant  strain  triangles  are  used  to  represent  a 
continuum,  and  the  hydrostatic  pressure  in  a  given  element  is  computed  using 
the  Mle-Grun-isen  equation  of  state.  The  code  treats  problems  involving  wave 
propagation  ind  elastic-plastic  flow.  It  is  arranged  to  provide  solutions  for 
projectile-target  impacts  and  explosive  detonation  problems. 

Material  properties  used  for  the  numerical  computations  are  given  in 
Table  1. 


Table  1.  Material  Properties 


Density  (kg/M^) 
Shear  Modulus  (GPa) 
Yield  Stress  (GPa) 
Ultimate  Stress  (GPa) 


Steel  (Mar-M-300) 

7855.00 

79-57 

2.07 

2.07 


Aluminum  (7Q75-T6) 

278*1 .00 
27.51 
0.29 
0.32 


The  propagation  of  elastic  wave3  in  a  3teel  bar  has  been  studied  to 
ascertain  t.ic  effects  of  mesh  orientation  and  artificial  viscosity  on  the 
numerical  solution  by  comparison  with  analytical  solutions.  The  analytical 
solutions  assume  a  3erai-infinite  bar.  The  calculations  presented  herein  are 
for  a  long  bap,  length-to-diameter  ratio  of  50.  Comparisons  between  the 
computational  results  and  the  analytical  solutions  are  made  before  the  arrival 
of  any  reflected  waves  from  the  rear  surface.  Thus,  the  finite  length  bar 
modeled  was  effectively  a  semi-infinite  bar  at  the  times  at  which  the 
comparisons  are  made.  The  elastic  impact  simulation  involves  a  bar  of 
circular  cross-section  made  of  high-hard  steel.  The  length  of  the  bar  Is  12.7 
cm  and  the  diameter  is  0.25*1  cm.  The  bar  impacts  a  rigid  wall  with  a  striking 
velocity  of  3.0*i8  m/sec. 

In  the  numerical  simulation  of  the  elastic  bar,  a  grid  of  quadrilateral 
elements  comprised  of  four  triangles  per  quadrilateral  w..s  used.  Various 
spatial  discretizations,  such  as  element  aspect  ratio  and  the  number  of 
elements  across  the  bar'3  radius,  have  been  studied.  The  elementary  one- 
dimensional  wave  equation  and  the  approximate  theory  for  determining  how 
elistic  3tr?ss  waves  propagate  in  infinite  bars,  given  by  Richard  Skalak,  are 
compared  with  the  computational  results  in  Section  III. 

The  elementary  one-dimensional  linear  wave  equation  Is 


2  2 
3  u  2  3  u 

2  “  c  2 
3t  32 


(1) 
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where  u  Is  the  displacement  along  the  •/.  axis  of  the  bar.  The  constant,  c,  is 
the  elementary  wave  velocity  for  a  bounded  medium  and  is  qiven  by 


and  Uiu  exact  dilatations!  wave  velocity  in  an  infinite  elastic  medium  is 


VEC1  • 
(1.- 


V)  (1  -  2V  )p  , 


where  »•  is  Poisson's  ratio,  E  13  Young's  Modulus,  and  /•  is  the  density. 
According  to  the  elementary  theory,  a  pulse  propagates  along  the  bar  without 
changing  its  shape.  The  general  form  of  the  wave  solution  is 


u  =  f(z  -  c,t)  +  g{z  +  c,t)  . 


(H) 


The  specific  shape  of  the  wave  is  determined  by  the  form  of  f.  Since  the  wave 
velocity,  c,  is  constant  for  an  ideal  onc-dimensional  elastic  medium,  the 
shape  of  the  wave  doc3  not  change  during  propagation.  The  equation  of  motion 
given  in  (1)  is  vaiid  for  wave  propagation  in  slender  bars  (length-to-dlameter 
ratio  of  at  least  10)  and  neglects  the  effects  of  transverse  strain,  lateral 
inertia,  body  forces,  and  dissipative  forces.  However,  the  equation  can  be 
applied  with  reasonable  success  to  problems  where  the  longest  wavelength  Is 
six  to  ten  times  greater  than  the  typical  cross-sectional  dimension  of  the 
bar. 


The  analysis  of  the  transmission  of  waves  in  a  uniform  circular  cylinder 
of  finite  length  i3  complex  because  of  the  Introduction  of  the  end  conditions 
on  the  bar.  Richard  Skalak  has  obtained  a  solution,  which  includes  transverse 
strain  effects,  for  elastic  wave  propagation  in  a  semi-infinite  elastic  bar 
subjected  to  a  suddenly  applied  impulse  through  the  superposition  of  two 
parts,  namely 

82U  fXy  +  2y)  34  3u>n 

— ~  - +  2y 

3t  p  3r  3x  ,  (5) 


32U  (X  +  2y)  34  2  y  3 

— 2“  - c™0) 

3t;  p  3z  r  3r  ’ 
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where 


1  £(rU  )  3Uz 

A  - - + - 

x  3r  3z 


(3u  3u  \ 

3z  3r  / 


J.  X  +  2y 
c  - - — 

P 


Skilak'3  30  Jtion  is  valid  for  late  times  at  points  far  away  from  the  point  of 
application  of  the  lond  and  13  obtained  in  terms  of  Airy's  integral.  The 
approximate  solution  is  of  the  form 

Ai(o)da  +  ^  Ai(a)daj  ,  (\o) 


where  »  is  hoisson  ratio,  «'  and  «"  are  the  non-dimensional  coordinates  which 
measure  th»*  distance  at  any  point  along  the  wavelength,  with  reipect  to 

velocity,  and  is  Airy '3  integral  with  respect  to«  ,  a  function  of 

time.  Fig  ire  l  depicts  Skalak's  refinements  to  the  elementary  theory 
solution.  Vhese  solutions  provide  a  guide  to  the  reliability  and  accuracy  of 
the  numeric  d  simulations. 

When  coaling  with  numerical  solutions  to  met  problems,  difficulties 
arise  when  >hock  waves  are  present.  Von  Neuman  and  Richtmyer  first  introduced 

artificial  viscosity*1  in  order  to  treat  problems  involving  physical  and 
mathematical  discontinuities  in  a  continuum.  Their  one  term  formulation, 


-  c2Q  p(Ax)2 


(lx) 


where  q1  =  0  for  dx/dx  greater  than  or  equal  to  zero  smears  the  shock  front 

over  several  mesh  widths,  thus  changing  the  discontinuities  into  steep  but 
continuous  wavefronts  which  can  be  treated  in  the  context  of  continuum 
mechanics. 

Later  on,  it  was  found  that  many  finite  difference  discretizations  tended 
to  introdur :  spurious  wave  motions  into  the  numerical  computations.  In  order 


Wilkins,  M.  L. ,  "Use  of  Artificial  Viscosity  in  Multi-Dimensional  Fluid  Dy¬ 
namic  Ca  eulations, "  Journal  of  Comp.  Physics ,  Vo l  36,  No.  6,  July  15,  1980. 


Figure  1.  Richard  Skalak's  Analytical  Solution. 


to  suppress  these,  the  artificial  viscosity  concept  of  Von  Neuman  and 
Ri^httnyer  was  extended  to  include  a  second  term, 


*  ■  “1  w 


a,  -  Ck/p(fc02  III  |  ('3> 

and  q^  r  0  for6x/dx  greater  than  or  equal  to  zero.  The  main  purpose  of  the 

second  term  was  to  suppress  numerical  noise  as  the  shock  wave  propagates.  For 
computations  involving  solids  where  the  sound  speed  exists  even  at  zero 
pressure,  the  linear  and  quadratic  terms  are  used  together.  Hence,  the 
generalized  fjrm  is 


q  »  kjpcAx  |||]  - 

where  c  is  the  sound  speed,  ax  Is  the  characteristic  dimension  in  the 
direction  of  wave  propagation,  pis  the  density,  k^  is  the  linear  component, 

k^  is  the  quadratic  component,  and  x  is  the  coordinate  in  the  direction  of 

motion.  In  EPIC-2,  the  generalized  formula  is  added  to  the  equation  of  state, 
which  helps  to  eliminate  numerical  overshooting  and  lower  level  noise  arising 
from  discontinuities. 

Having  established  the  effects  of  irregular  meshing  and  artificial 
viscosity  on  problems  of  elastic  wave  propagation,  an  additional  set  of 
simulations  was  undertaken  to  obtain  an  appreciation  of  the  effect  of  these 
parameters  on  situations  involving  plastic  flow.  Only  the  mesh  orientation 
was  varied  in  the  plastic  impact  study.  Three  mesh  orientations  were 
considered  which  are  shown  in  Figure  2. 

The  effects  of  the  various  mesh  orientations  can  be  seen  from  the 
deformation  patterns  but  become  even  more  evident  when  strain  and  pressure 
distribution::  are  studied.  Finally,  while  not  explicitly  considered  in  this 
study,  the  behavior  of  numerical  results  in  the  presence  of  different  failure 
criteria  and  post-failure  models  can  be  inferred  from  the  results  presented. 


3x 

3x 

3x 

3x 

(Ill) 


III.  RESULTS 


A.  Elastic  Impact 

A  comparison  of  the  results  from  the  various  computations  to  determine 
the  grid  size,  using  quadrilateral  elements  comprised  of  four  triangles,  gave 


IDIA  =  1  Hypotenuse  of  the  Triangle  is  drawn  from  the 
lower  left  hand  corner  to  the  upper  right  ha 
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Fin^ire  2.  Elenent  Orientation::  in  F.F1C-2  Co*- 


some  Interesting  results,  as  shown  in  Figures  3  -  l»l .  In  all  six  cases,  the 
change  in  amplitude  and  time  remained  relatively  small.  Using  a  uniform  grid 
and  refining  the  mesh,  there  is  one  considerable  effect:  by  increasing  the 
number  of  nodal  points  and  decreasing  the  size  of  the  elements,  the  higher 
frequencies  are  resolved. 

In  Figures  *l ,  6  and  8  the  wave  profile  in  the  bar  at  20ms  i3  presented 
for  l,  2  and  3  elements  across  the  radius  respectively.  In  all  three  cases 
the  olcnent  aspect  ratio  was  a  constant  of  1:1.  Thus,  as  the  number  of 
elements  across  the  radius  is  increased,  the  resolution  is  improved  since  the 
higher  frequency  waves  are  now  being  accounted  for  in  the  numerical  solution. 
Comparison  o°  the  high  resolution  wave  profile  in  Figure  8  with  Skalak's 
solution  in  Figure  l  shows  close  agreement. 

In  cont -ast  to  the  elementary  theory  solution,  the  wavefront  predicted 
from  both  Skalak's  solution  and  the  numerical  solution  has  a  finite  3lopc. 
The  elementary  theory  reasonably  approximates  some  features,  such  as  peak 
amplitude  of  the  behavior  of  long  bars  subjected  to  impulsive  loading;  but  it 
is  Incapable  of  resolving  some  aspects  of  the  problem.  The  solution  implies 
that  the  waves  transmitted  are  traveling  with  constant  velocity  and  without 
alteration  of  shape.  In  the  Skalak  and  numerical  solutions,  the  oscillations 
arise  due  to  transverse  strain  effects.  They  are  not  artificially  introduced 
by  the  numerical  scheme.  Furthermore,  as  one  refines  the  mesh,  the  wave 
profile  exhibits  on  increase  in  the  number  of  oscillations  behind  the  wave 
front. 


Figures  8,  10,  12  and  l»)  show  the  wave  profile  in  the  bar  at  20ms  for 
element  aspect  ratios  of  1:1,  2:1,  3:1  and  H:1  respectively.  As  the  aspect 
ratio  is  increased  the  axial  resolution  decreases.  Thi3  decrease  in  axial 
resolution  prevents  the  numerical  solution  from  capturing  the  high  frequency 
wave  components  in  the  axial  direction.  This  can  be  observed  in  the  plots  by 
the  decrease  in  the  number  of  oscillations  in  the  wave  profile  behind  the  wave 
front. 

The  cost  of  computations  may  often  be  a  major  factor  in  determining  the 
grid  size.  With  a  large  aspect  ratio  and  a  coarse  grid,  the  CPU  time  and 
costs  are  relatively  low.  While  the  cost  and  the  CPU  time  practically  double 
with  grid  refinements.  Thus,  the  user  must  make  a  trade  off  between  accuracy 
and  cost. 

The  effects  of  changes  in  artificial  viscosity  coefficients  are  shown  in 
Figures  15  -  17.  An  increase  in  the  linear  component  results  in  excessive 
distortion  of  the  wave  front  end  a  decrease  in  the  oscillations  occurring 
along  the  wavelength.  The  effect  of  changes  in  the  quadratic  coefficient  is 
small.  This  is  primarily  because  the  quadratic  component  dissipates  more 
rapidly  away  from  the  shock  than  the  linear  component. 


B.  Plastic  Impact 

An  additional  set  of  simulations  has  been  performed  to  gain  an 
appreciation  for  the  effect  of  element  orientation  on  the  numerical  solution 
of  impacts  Involving  plastic  flow.  In  this  set  of  computations  the  mesh 
orientation,  in  the  target  only,  was  altered.  All  computations  were 
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Figure  3.  Computational  Grid  Map,  '100  Elements. 
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terminated  at  lO/»s,  at  which  time  the  problems  are  sufficiently  advanced  to 
obtain  reasonable  results.  The  plastic  impact  studied  consists  of  a  steel 
sphere  impacting  an  aluminum  plate  with  a  velocity  of  152*1  m/sec.  The  target 
plate  thickness  is  5.08  cm.  The  first  target  mesh  orientation  studied  ha3  the 
diagonals  oriented  with  a  positive  slope  (IDIA=1  in  Figure  2),  see  Figure  18. 

Figures  19  and  20  show  the  computational  results  at  5;»3  and  10;is 
respectively.  At  10/is  the  target  is  forming  a  V-ahape  which  will  produce 
unrealistic  deformations  of  the  sphere  and  numerical  Instabilities.  The 
computational  results  for  the  diagonals  oriented  with  a  negative  slope  are 
presented  in  Figures  21  and  22.  At  5 /‘a  there  is  a  slight  pulling  of  the 
elements  in  a  positive  direction.  By  IO/13  the  target  deformation  is 
unreasonably  stiff.  Clearly,  from  Figures  19-22,  the  diagonal  orientation  ha3 
a  significant,  influence  on  the  predicted  deformations. 

In  Figures  23  and  2*1  the  quadrilateral  elements  comprised  of  four 
triangles  allow  for  largo  deformations  to  occur  without  introducing  artificial 
stiffness  and/or  asymmetries  into  the  calculations,  which  are  very  typical  of 
the  triangular  zones. 

The  response  of  the  target  material  has  been  carefully  examined  for 
strain,  pressure,  and  distribution  as  well  as  failure  patterns.  Material 
failure  in  the  EPIC-2  code  is  governed  by  a  failure  criterion  based  on 
equivalent  plastic  strain.  When  a  user’s  specified  value  for  this  quantity  in 
.  an  element  is  exceeded,  the  element  is  considered  to  have  partially  failed, 
and  no  tensile  or  shear  stresses  are  permitted.  In  effect,  the  element 
behaves  as  a  fluid.  Total  failure  at  higher  values  of  equivalent  plastic 
strain  is  also  permitted.  In  this  case,  neither  stress  nor  pressure  can  be 
carried  by  an  element,  and  it  is  therefore  ignored  in  the  computational  loop. 
The  element  connectivity  is  destroyed>  but  the  nodes  and  their  associated  mass 
and  velocity  are  retained  to  conserve  mass  and  momentum.  The  results 
pertaining  to  the  failure  of  the  elements  in  the  target  are  shown  in  Appendix 
A.  Those  elements  marked  with  an  "F"  have  exceeded  the  partial  failure 
criterion  of  0.6. 

With  regard  to  the  first  two  cases,  107  (IDIA=l)  and  108  (IDIA=2),  the 
elements  are  being  failed  far  away  from  the  impact  area.  The  total  amount  of 
failure  is  small;  and  the  dispersement  of  the  failed  elements  is  more  widely 
distributed  in  the  target,  in  comparison  to  Case  109* 

Element  strains  are  derlyed  by  integrating  the  strain  rates,  which  are 
obtained  from  the  geometry  of  the  elements  and  the  nodal  velocities.  The 
strain  contours  are  evaluated  at  five  different  levels  and  are  shown  in 
Appendix  B.  The  strain  contours  in  Cases  107  and  108  sharply  deviate  from  the 
area  of  impact  and  extend  out  farther  into  the  target  in  comparison  to  Case 
109. 


The  hydrostatic  pressure  is  dependent  on  the  volumetric  strains  and 
internal  energy  in  the  elements.  The  results  from  the  positive  and  negative 
pressure  cjntours  are  shown  in  Appendix  C.  In  comparing  the  first  two  Cases 
to  the  third,  the  positive  pressure  contours  are  more  condensed  around  the 
area  of  impact  and  the  failed  elements.  With  regard  to  the  negative  pressure 
fields,  there  is  a  greater  magnitude  of  the  pressure  in  the  projectile  than  in 
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Figure  23.  Deformation  of  Sphere/Plate  Impact,  t 
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the  target,  and  the  majority  or  the  pressure  fields  in  the  target  are  located 
in  very  small  regions. 


IV.  CONCLUSIONS 


The  elastic  impact  analyses  demonstrate  that  close  agreement  between  the 
analytical  (Skalak)  and  numerical  solutions  can  be  obtained  provided  a  fine 
grid  with  an  aspect  ratio  of  1:1  is  used.  Coarse  computational  grids  and  high 
aspect  ratios  prevent  the  numerical  solution  from  capturing  the  high  frequency 
components. 

For  the  effects  of  artificial  viscosity  upon  wave  propagation  problems, 
the  linear  term  provides  the  additional  dissipation  required  to  obtain  a 
smooth,  continuous  solution,  but  it  also  tends  to  broaden  the  shock  front, 
the  linear  component  should  range  between  the  values  of  0.05  to  0.20  for  the 
oest  results.  In  most  computer  programs,  the  linear  and  quadratic  components 
are  used  in  conjunction  with  each  other  in  order  to  stabilize  the  effects  of 
excessive  numerical  ringing.  The  quadratic  component  is  most  effective  when 
the  values  range  between  2.0  and  M .0  for  the  purpose  of  obtaining  accurate 
results. 

In  the  cases  studied,  for  the  high  velocity  impact  of  a  sphere  into  a 
flat  plate,  the  effects  upon  the  failure  of  the  elements,  the  deformation 
behavior  characteristics  of  the  elements,  and  the  pressure  and  strain  fields 
show  that  quadrilateral  elements  comprised  of  four  triangles  per  quadrilateral 
Is  the  most  reasonable  representation  of  observed  behavior  in  elastic-plastic 
impact  problems. 
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Figure  C5.  Pressure  Contour  top,  t  =  10/is,  Case  10B 


Figure  C 6.  Pressure  Contour  Map,  t  =  lO^s,  Case  109 
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Wilmington,  MA  01887 

3  BntccLlc  Columbus  Laboratories 
505  King  Wenuc 
ATTN:  Dr.  G.  T.  Hahn 

Dr.  L.  E.  Hulbect 
Dr.  S.  Snmpnth 
Columbus,  OH  4320 1 

A  Boeing  C  mpany /Aerospace  Division 
ATTN:  Mr.  R.  G.  Blaisdell 
( M . S .  A0-25) 

Dr.  N.  A.  Armstrong, 

C.  J.  Artura 
(M.S.  8C-23) 

Dr.  B.  J.  Henderson 
(M.S.  43-12) 

P.O.  Box  3707 
Seattle,  WA  9G12A 

2  Brunsulck  Corporation 
A300  Industrial  Avenue 
ATTN:  P.  S.  Chang 

R.  Grover 
Lincoln,  ME  6850A 

1  Computer  Code  Con¬ 
sultants,  Inc. 

1680  Camlno  Redondo 
ATTN:  Dr.  Wally  Johnson 
Los  Alamos,  NM  875AA 

1  Dresser  Center 
PO  Box  1A07 

ATTN:  Dr.  M.  S.  Chawla 
Houston,  TX  77001 

1  Effects  Technology,  Inc. 

5383  Hollister  Avenue 
Santa  Barbara,  CA  93111 


1  Electric  Power  Research 
Institute 
PO  Box  10412 

ATTN:  Dr.  George  SI iter 
Palo  Alto,  CA  94303 

l  FMC  Corporation 

Ordnance  Engineering  Division 
San  Jose,  CA  951 1A 

1  Ford  Aerospace  and 

Communications  Cor¬ 
poration 

Ford  Road,  PO  Box  A 
ATTN:  L.  K.  Goodwin 
Newport  Beach,  CA  92663 

l  General  Dynamics 
PO  Box  2507 
ATTN:  J.  H.  Cuadros 
Pomona,  CA  91766 

1  General  Electric  Company 
Lakeside  Avenue 
ATTN:  D.  Graham,  Rm  1311 
Burlington,  VT  05A01 

1  President 

General  Research  Corporation 
ATTN :  Lib 
McLean,  VA  22101 

1  Goodyear  Aerospace  Corporation 
1210  Massilon  Road 
Akron,  OH  AA315 

1  H.  P.  White  Laboratory 
31 l A  Scarboro  Road 
Street,  MD  2115A 

5  Honeywell,  Inc. 

Government  and  Aerospace 
Products  Division 
ATTN:  Mr.  J.  Blackburn 
Dr.  G.  Johnson 
Mr.  R.  Simpson 
Mr.  K.  H .  Doeringsfeld 
Dr.  D.  Vavrick 
600  Second  Street,  NE 
Hopkins,  MN  553A3 
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Hughes  Aircraft  Corporation  1 

ATTN:  Mr.  W.  Keppcl 

MS  M-5,  Bldg  808 
Tucson,  A7.  85706 

Kaman  Sciences  Corporation  2 

1500  Garden  of  the  Gods  Road 
ATTN:  Dr.  P.  Snow 

Dr.  D.  Williams 
Colorado  Springs,  CO  80907 

Lockheed  Palo  Alto  Research  1 

Laboratory 
3251  Hanover  Street 
ATTN:  Org  5230,  Bldg  201 
Mr.  R.  Robertson 

Palo  Alto,  CA  94394  * 

Lockheed  Missiles  and  Space 
Company 
PO  Box  504 

ATTN:  R.  L.  Williams 

Dept.  81-11,  Bldg  154  : 

Sunnyvale,  CA  94086 

Materials  Research 
Laboratory,  Inc. 

I  Science  Road 
Glenwood,  IL  60427 

McDonnell-Douglas  Astro¬ 
nautics  Company 
5301  Bolsa  Avenue 
ATTN:  Dr.  L.  B.  Greszczuk 
Dr.  J.  Wall 

Huntington  Beach,  CA  92647 

New  Mexico  Institute  of 
Mining  and  Technology 
ATTN:  TERA  Group 
Socorro,  NM  87801 

Northrop  Norair 
3901  W.  Broadway 
ATTN:  R.  L.  Ramkumar 
Hawthorne,  CA  90250 


Organization 

Nuclear  Assurance  Corporation 
24  Executive  Park  West 
ATTN:  T.  C.  Thompson 
Atlanta,  GA  30245 

Orlando  Technology,  Inc. 

P0  Box  855 

ATTN:  Mr.  J.  Osborn 
Mr.  D.  Matuska 
Shnlimar,  FL  32579 

Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSQ-OR 
Fort  Bcnning,  GA  31905 

Rockwall  International 
Missile  Systems  Division 
ATTN:  A.  R.  Glaser 
4300  E.  Fifth  Avenue 
Columbus,  OH  43216 

Schumberger  Well  Services 
Perforating  Center 
ATTN:  J.  E.  Brooks 
J .  Brookman 
Dr.  C.  Aseltine 
PO  Box  A 

Rosharon,  TX  77543 

Science  Applications,  Inc. 

101  Continental  Boulevard 
Suite  310 

El  Segundo,  CA  90245 
Commander 

US  Army  Development  &  Employment 
Agency 

ATTN:  MODE-TED-SAB 
Fort  Lewis,  WA  98433 

S-CUBED 
PO  Box  1620 

ATTN:  Dr.  R.  Sedgwick 
La  Jolla,  CA  92038-1620 
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2  TRW 

One  Space  Park,  Rl/2120 
ATTN:  D.  A'jshcrman 
M.  Bronstcin 
Redondo  Beach,  CA  90278 

1  United  Technologies 
Research  Center 
438  Weir  Street 
ATTN:  P.  R.  Fitzpatrick 
Glastonbury,  CT  06033 

1  US  Steel  Corporation 
Research  Center 
125  Jamison  Center 
Monroeville,  PA  15146 

1  VPI  &  SU 

106C  Norris  Hall 
ATTN:  Dr.  M.  P.  Knnat 
Blacksburg,  VA  24061 

2  Vought  Corporation 
PO  Box  225907 
ATTN:  Dr.  G.  Hough 

Dr.  Paul  M.  Kenner 
Dallas,  TX  75265 

1  Westinghouse,  Inc. 

P0  Box  79 

ATTN:  J.  V.  Fan 

W.  Mifflin,  PA  15122 

1  Drexel  University 

Department  of  Mechanical  F.ngr 
ATTN:  Dr.  F.  C.  Chou 
32d  and  Chestnut  Streets 
Philadelphia,  PA  19104 

Southwest  Research  Institute 
Dept,  of  Mechanical  Sciences 
ATTN:  Dr.  U.  Lindholm 
Dr.  W.  Baker 
Dr.  R.  White 
Dr.  M.  F.  Kanninen 
Dr.  C.  Anderson 
8500  Culebra  Road 
San  Antonio,  TX  78228 


4  SRI  International 

333  Ravenswood  Avenue 
ATTN:  Dr.  L.  Seaman 
Dr.  L.  Curran 
Dr.  D.  Shockey 
Dr.  A.  L.  Florence 
Menlo  Park,  CA  94025 

2  University  of  Arizona 

Civil  Engineering  Department 
ATTN:  Dr.  D.  A.  DaDeppo 
Dr.  R.  Richard 
vcson,  A 7.  85721 

„  versity  of  Arizona 

.iLool  of  Engineering 
ATTN:  Dean  R.  Gallagher 
Tucson,  A 7.  85721 

1  University  of  California 
Los  Angeles 
ATTN:  Dr.  N.  7.iv 
Los  Angeles,  CA  90024 

1  University  of  California 
Department  of  Physics 
ATTN:  Dr.  Harold  Lewis 
Santa  Barbara,  CA  93106 

2  University  of  California 
College  of  Engineering 
ATTN:  Prof.  W.  Goldsmith 

Dr.  A.  G.  Evans 
Berkeley,  CA  94720 

2  University  of  Delaware 
Department  of  Mechanical 
Engineering 

ATTN:  Prof.  J.  Vinson 
Prof.  B.  Pipes 
Newark,  DE  19711 

1  University  of  Denver 

Denver  Research  Institute 
ATTN:  Mr.  R.  F.  Recht 
23 SO  S.  University  Blvd. 
Denver,  CO  80210 
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2  University  of  Florida 

Department  of  F.ngincering 
Sciences 

ATTN:  Dr.  P,.  L.  Sicrakowski 
Dr.  L.  F..  Malvern 
Gainesville,  FL  .52601 

l  University  of  Oklahoma 
School  of  Aerospace, 
Mechanical  and  Nuclear 
Engineering 
ATTN:  Dr.  C.  W.  Bert 
Norman,  OK  73069 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
\TTN :  A-7SY-D 

A  iXSY-MP ,  II.  Cohen 
(Mr,  USATECOM 
ATTN:  AMSTK-TO-F 
l..lr,  USACSTA 
\TTN :  Mr.  W.  Plcss 

Mr.  S.  Keithley 

Cdr,  CROC,  AMCCOM 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  commcnts/answcrs  to  the  items/questions  below  will 
aid  us  in  our  efforts. 


1.  BRL  Report  Number _ 

2.  Date  Report  Rcccived_ 


Date  of  Report 


3.  Docs  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 
other  area  of  interest  for  which  the  report  will  be  used.) _ 


<1.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


OLD 

ADDRESS 


Organization 


Address 


City,  State,  Zip 


(Remove  this  sheet  along  the  perforation,  fold  as  indicated,  staple  or  tape 
closed,  and  mail.) 


•>  v  rr  •ej <v  ■ 


FOLD  HERE 


Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-0D-5T 

Aberdeen  Froving  Ground,  MD  21005-5066 
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Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-OD-ST 

Aberdeen  Proving  Ground,  MD  21005-9989 


FOLD  HERE 


